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abstract 


SeVtral meth ° ds ° f transition prediction by U„ear 
stability analysis are compared. The spectral stability 
analysis code SALLY is used to analyse now* over laminar 
now control wi„ 9s . it is shown that transition prediction 
by the envelope r-thod and a new modeled wave packet method 

are comparable in reliability but that th« i 

y uc that the envelope method 

is more efficient computationally. 
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maximum disturbance amplitude 
Chebyshev coefficients 
wing chord 

dimensional frequency 
wave number vector 
algebraic mapping parameter 
N-factor « ?n A/A Q 

displacement thickness Reynolds number, U^6 Vv^ 

chord Reynolds number, c/v^ 

arc length aiong an arbitrary path on the wing 

Chebyshev polynomial 

time 

unperturbed x-velocity in the boundary layer 
potential flow vector at edge of boundary layer 
x-component of 

incoming free stream velocity 
group velocity vector 

perturbation velocity in the y-direction 
unperturbed 2-velocity in the boundary layer 
mapped coordinate normal to wina surface 
coordinate in the direction of the normal chord 
coordinate normal to the wing surface 
coordinate along the wing span 


li 


a x-wave number 

a ft angle of attack 

8 2 -wave number 

oi frequency (complex) 

ft frequency (real) 

* 

f> displacement thickness 

X wave length 

v kinematic viscosity 

0 wing sweep angle 

ip angle formed by the wave number vector with the x-axis 
^g angle formed by the group velocity vector with the x-axis 

'I'p angle formed by the p itential flow vector with the x-axis 

$ eigenfunction; defined in Eq. (3) . 

n adjoint of the e 4 -»n function 
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SECTION 1 - INTRODUCTION 


In this report several methods of transition prediction 
using linear stability analysis are compared. The incompress 
lble linear stability computer code SALLY is used in various 
ways to study three-dimensional boundary layer flow over 
laminar flow control (LFC) wings. Here we compare the so 
called envelope method 1 with wave-packet methods 2 , to predict 
transition. We conclude that the envelope method is at least 
as reliable as the more complicated and less efficient wave 


packet method. 

Consider the stability of three dimensional laminar flow 
over swept wings with sweep angle 0. The coordinate system 
used on the wing is depicted in Fig. 1. The x-axis is in the 
direction of the normal chord, the y-axis is normal to the 
surface of the wing while the z-axis is along its span. 

Neglecting the curvature of the wing surface, compressi- 
bility effects, and non-parallel flow effects, linear distur- 
bances satisfy the Orr-Sommer feld equation 


- « 2 - bV* 
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iR { (aU + BW - in) [ — r 
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with the boundary conditions 


*(0) * (0) * 0 ; <M°°) bounded. (2) 

Here the perturbation velocity in the y-direction is assumed 
to be of the form 

. _ fx/ . i (ax+Bz-uit) , 

v' = Re[$(y)e 1# (3) 

0(y) and W(y) are the (unperturbed) laminar boundary 
layer velocities in the x- and z-directions, respectively, 
and R is the Reynolds number. It xS assumed that all 
variables are non-dimensionalized with boundary layer 
scaling. 

Equations (1) - (3) constitute an eigenvalue 
problem for the frequency w and wavenumbers ct,3 . 

For given Reynolds number R, this eigenvalue problem 
provides a complex dispersion relation of the form 

ts) = a>(a,6) (4) 

relating the complex parameters a,B and w. 


... # 


1*2 





Semi-empirical methods to predict transition on 

LFC wings are based on tracing the evolution of modes 
1 

across the wing. An appropriate N— factor for transition 
correlation is defined as the (logarithm of the) total 
growth factor across the wing (see below). A good transition 
predictor is one for which transition occurs at nearly constant 
N for a wide variety of wings and flow conditions. 

For natural transition, disturbances of all frequencies 
are present on the wing surface. In this case, there are 
many optional ways to compute N factors. The first choice 
is between temporal and spatial stability theory. In 
temporal theory, a and 6 are real while w is complex? 
the mode grows in time if Im(w) > 0, but the mode does 
not grow in space. An N- factor for transition correlation 
may be defined as 

N ■ | Im(u>)/|Re (v g ) |ds 

where v = (3o)/3 a, 3u>/38) is the (complex) group velocity 

g 

and s is the arclength along an appropriate curve on the wing. 
The N-f actor (5) is not fully defined until a prescription 
is given for singling out a specific mode at each position 
on the wing and for defining a specific curve on which to 
integrate. We shall return to these questions in Sec. 2. 

In spatial stability theory, w is real but a and/or 
0 may be complex. Again, there is arbitrariness in the 
definition of an appropriate N- factor because of the variety 
of excitable modes on the wing. 



SECTION 2 


WAVE PROPAGATION IN BOUNDARY LAYERS 

The complex eigenvalue relation (4) provides 

two real relations among the three complex quantities 

o,0, and u>. In temporal stability theory, the 

requirements that a and 6 be real provide two more 

conditions so there remain two arbitrary parameters 

among Re (a) , Re (6) , Re(u») , and Im(<i>) . 

There are several ways to remove this arbitrariness 

in the computation of the growth factors N. In the 

envelope method 1 , Im(uj) is maximized with respect 

to a at fixed Re u> Iwhich then determines o,0 

and u) uniquely at each point on the wing] and the 

curve in (5) is defined to be everywhere tangent 

to Re (v) . 

9 

With spatial stability theory, there remain 

three independent real parameters among a, 6 and 

Re(iu) once the eigenvalue condition (4) is satisfied. 

One possibility is to require that the direction of 

most rapid growth, which is parallel to the vector 

(-lm(a), -Im (0) ) , be parallel to Re(v g ) and that 

the resulting value of the most rapid growth rate 

be maximized with respect to the remaining two 

3 

independent parameters. 

Alternatively, it is possible to use wave 
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packet theory to remove the arbitrariness in the definition 
of N- factors. For a conservative dynamical system, 
kinematic wave theory implies that a wave packet 

propagates in physical and wavevector space according 

4 

to the Hamilton- Jacobi equations 


dx 

- 3(0 

(6) 

dt. 

da 

dz 

3a> 

(7) 

dt 

33 

da 

_ _ 3u> 

(8) 

dt 

3x 


_ _ 3ii) 

(S) 

dt 

3z 


Nayfeh considered the extension of Eqs. (6)- 
(9) to non-conservative systems where a,$, and w can 
Jae complex. Then, and may also be complex. 

For a physical solution with real x,z, and t to 
exist, (6) and (7) imply that the group velocity 
(u> ,(d.) must be real. Nayfeh proposed the computation 

Ot p 

of '.wave packet solutions determined by the six independent 

conditions: (i) the eigenvalue condition (4); 

(ii) Im u) = Im (d 0 = 0; (iii) Re fixed; (iv) 
a 0 

Re 6 $ixed; and (v) dx/dt = dz/dt = Under 

these conditions the N-factor is determined by 
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t 

N * / I-w a Im(a)-Wglm(6)-»-ImCu) Jdt (10) 

t 0 

Finally we study a modified non-conservative wave 

packet formulation in which a, 6* and o> are 

determined by: (i) the eigenvalue condition (4); 

(ii) Im (d = Im <u . = 0; (iii) Rem fixed with Im u» = 0; 
a p 

and (iv) dx/dt = oj , d z/dt = w. . The motivation for 

a d 

these latter conditions is simply that laminar flow 
over a LFC wing may be assumed steady so a wave packet 
should propagate at fixed real frequency. The N- 
faetor is given by (10) with Im(u>) = 0. 

Calculations made with Nayfeh's formulation of 
the wave packet equations were extremely sensitive 
and gave transition prediction at highly variable end 
unpredictable values of the N factor. Therefore, we 
report in Sec. 4 only the results obtained by our 
modified wave packet formulation in which the condition 
Re 0 fixed is dropped in favor of Im<i>= 0. 


SECTION 3 - NUMERICAL METHOD 


In the computer code SALLY 1 , Eqs. (1) - (2) arc solved 
using a spectral method based on Chebyshev polynomials 5 . The 
boundary layer direction y 0 <_ y < «>, i s mapped into the finite 
interval -1 w < 1 by the algebraic mapping 


w 

and 4>(y) 
series 


= 2 
is 


- 1 


y+L 
approximated 


(ID 

as the finite Chebyshev polynomial 


M 

<Hy) = E a n T (w) (12) 

n=0 


The resulting algebraic eigenvalue problem is solved globally 
(if a guess for the eigenvalue is not available) by a generalized 
QR algorithm or locally (if a good guess is available) by inverse 
Rayleigh iteration**. The resulting scheme is very efficient 
and accurate. 

^ i 

Group velocity (Vg = — — ) can be calculated using the adjoint 

3k 

eigenfunction oi the Orr-Sommerfeid equation. Thus if the Orr- 
Somnrierfeld equation is written as 

L (k , x, u) (k) ) <J> = 0 (13) 


then the differentiation with respect to wave number k gives 


3L . 

— ♦ + 

3L 

^ * * 

t* 

j V 
4-l-e- 
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3k 
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3k 

3k 


(14) 
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Taking the inner product of (14) with the adjoint n of 


the eigenfunction <J> gives 


3d) 

9k 


(n,% <J>) 

9k 

/ 31 *.* 

(n ' 35 * > 


(15) 


In order to correlate transition we calculate N-factor for 
a fixed real frequency and repeat the calculations for several 
other frequencies to find the most dangerous mode. The eigen- 
value solver in SALLY 1 searches for the eigenvalue w given a 
and 8- So we iterate on Re a, Im a, Re 8 and Im 8 such that 
the conditions described above for modified wave packet method 
are satisfied. This is achieved usually in three iterations. 
Following equations are employed in this iterative process. 
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The properties of the laminar boundary layer profiles 
required to solve (1) — (2) are obtained using a compressible 
boundary layer code for swept tapered wings developed by Kaups 
and Cebeci 1 . 
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SECTION 4- 
RESULTS 


7 

Burrows has reported flight transition data 
taken at Cranfield for a large, untapered, 45° swept 
half wing mounted as a dorsal fin upon the mid-upper 
fuselage of an Avro Lancaster airplane. The airfoil 
section was made-up of two semi-ellipses,, one of which 
constituted a faired trailing edge and the oth r 
corresponding to the leading edge of a 10 percent 
thick airfoil, with effective chord of 10.83 feet, 
measured in the free stream direction. The location 
of the beginning of transition in the Cranfield data 
was estimated as given in Ref. 8 - Two of the Cranfield 
flight tests were chosen for correlating transition 
using wave packet theory* 

In the first test, case, calculations were made 
for a chord Reynolds number of 11 . 7 x 10 ® and -2° 
angle of attack. In this flow, transition begins c*t 
x/c - 5.5%. A maximum N factor of 7.6 was obtained 
at a frequency of 1250 Hz both with the envelope 
method and the modified wave packet method, 

-the predicted variation of the N factor up 
to the transition location was almost identical for the 
envelope method and the modified wave packet method. 

We also compute the solution of the conservative wave 
packet equations (6 )-(*>) tn which only the real 
parts of equations (6) - (/) are taken while (8) -(9) 
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are solved in their full complex form. The resulting 
N factor at transition is 5.2. The variation of N 
factor with x/c for the various methods is plotted 
in Fig. 2. 

Wave angle, wave length and the direction of 
the group velocity as predicted by the envelope and 
wave packet methods are given in Figs. 3-5. Although 
the results are qualitatively similar, there is 
appreciable quantitative difference in these parameters 
at the transition location. It is surprising that the 
N factor calculated by the envelope and modified wave 
packet methods are the same. 

In the second test case, the angle of attack 
of the wing was changed to zero. In this case, 
transition occurred experimentally at x/c ~ 7%. 

The envelope method gave an N factor of 10.8 at 
a frequency of 1000 Rz. The wave packet method gave 
a maximum N factor of 10.5 at a frequency of 
1200 Hz, which is close to the prediction of the 
envelope method. The variation of N factor with 
x/c is plotted in Fig. 6. The predictions of the 
conservative wave packet approximation and a fixed 
wavelength, fixed frequency integration are also 
plotted in this figure. The conservative wave packet 
approximation gave an N factor at transition of 
8.6 rather than 10.5. 


Figure 7 shows the influence of frequency 
on N factor at transition for the wing as 
predicted by the wave packet theory. Wave angle, 
wave length and direction of the group velocity 
for this particular wing are shown in Figs. 8-10. 
Again there is substantial quantitative difference in 
the predictions of the two methods. 


SECTION 5 ~ CONCLUSIONS 


Calculations were made for a Cranfield 45° swept wing 

with Re = 11.7 x 10 6 using a modified wave packet method and 
c 

the envelope method. Both methods gave an N factor of 7.6 at 
transition location for an angle of attack, a p i ~ ~ 2°. For 
a = 0°, the envelope and modified wave packet methods gave N 

A 

factors of 10.8 and 10.5, respectively. Since it may be argued 
that the wave packet method is physically more relevant for pre- 
dicting transition in three dimensional boundary layers, it was 
initially hoped that the wave packet method might give more con- 
sistent transition N factors. However, the results show that 
the wave packet method provides N factors which are at best as 
consistent as those of envelope method. Since the wave packet 
method is at least 3 times as expensive to use as the envelope 
method, the latter is recommended for engineering design cal- 
culations . 

We would like to thank D. M. Bushnell and J. N. Hefner for 
helpful discussions. This work was supported by the National 
Aeronautics and Space Administration under Contract NAS1-15604 
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FIGURE CAPTIONS 


Figure 1. A plot of the coordinate system on a swept 
wing. 

Figure 2. A plot of N versus percent of chord 
x/c for various methods applied to a swept wing 
of an Avro Lancaster airplane at -2° angle of 
attack. Solid curve: modified wave packet method 

and envelope method at f = 1250 Hz which gives nearly 
the ma ximum N at the transition point. Dashed 
curve: result of integrating equations (6) - (9) 

across the wing with (6) and (7) replaced by 
their real parts. The curves are plotted from the 
beginning of the unstable flow region until the 
transition point at x/c = 5.5%. 

Figure 3. A plot of wave propagation angle versus 
x/c for the same flow as in Figure 2. 

Figure 4. A plot of wavelength versus x/c for the 
same flow as in Figure 2. 

Figure 5. A plot of the direction of the group velocity 
for the same flow as in Figure 2. 

Figure 6. Same as Fig. 2 except for the wing at 0° 
angle of attack. In addition to the results of the 
wave packet methods and envelope method t the N 
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factor obtained by integrating a fixed 
wavelength, fixed frequency mode across the wing is 
given. Here N is given by (5) and the mode is 
determined by the six real conditions: (i) F i4 

(ii) Im a = Im B = 0; (iii) X/c = 0.001; ( , 

Re a) = 7 50 Hz . 

Figure 7. Variation of N at transition versus 
frequency obtained using the modified wave 
packet method for the same flow as in Figure 6. 

Figure 8. A plot of wave propagation angle versus 
x/c for the same flow as in Figure 6. 

Figure 9. A plot of wavelength versus x/c for the 
same flow as in Figure 6. 

Figure 10. A plot of the direction of the group 
velocity for the same flow as in Figure 6. 
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